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Abstract 

We construct the list of all possible inert states of spin Bose condensates for S < 4. In doing so, 
we also obtain their symmetry properties. These results are applied to classify line defects of these 
spin condensates in zero magnetic field. 
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I. INTRODUCTION 



Optical trapping of atomic Bose gases has provided us with excitingnew 



physica 
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systems, which have attracted significant theoretical interests (e.g. 



xcitmgnew 



In constrast to 4 He, these Bosonic atoms possess finite hyperfine spins 
(hereafter simply as spins, denoted by S), hence internal degrees of freedom. Therefore in 
the Bose condensed phase, in contrast to 4 He where only the gauge symmetry is broken, 
the state can also have broken spin-rotational symmetries (we ignore the case of fragmented 
states l^] here). The present paper is devoted to symmetry considerations of these spin- 
condensates. Knowing the broken symmetries, as well as the residual symmetries, is essential 
in understanding the properties of these spin superfluids, in particular topological excitations 
such as vortices. 

The states that are expected to be realized are determined by minimization of (free) 



energy of the system. Such a program has been carried out by numerous authors 
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,13, 



131 ]. assumingthat mean-field theory applies. The minimization of this free energy is 



rather simple 



3| for S 



indeed rather complicated [12 



ready not quite trivial [8j] (see also [l5]) for S = 2, and 
13j for S = 3. There are also a few attempts to find out the 
residual symmetries of these minimum energy states found 16|, Il7 |. 

On the other hand, symmetry considerations can actually help in finding a special set of 
possible states known as "inert" states for the systems. These states are those where the 
wavefunction characterizing the state is actually independent of the (interaction) parameters 
that enter in the free energy that is being minimized. The form of these wavefunctions follow 
from symmetry considerations alone, depending only on the symmetries of the uncondensed 
(normal) state (thus independent of whether mean- field theory has been assumed). ^From 
such considerations, the symmetries of these inert states are thus also directly determined. It 
is however, true that this procedure only determines the possible stationary points (relative 
minima, extrema or saddle points) of the free energy functional. Additional calculations (in 
particular the actual form of the free energy and the parameters which enter it) are needed 
to determine which state is actually a stable minima. 

Such a program has actually been carried out before in other systems. In superfluid 3 He, 
Cooper pairing occurs in angular momentum L = 1 and spin S = 1 channel. Spatial and/or 
spin rotational symmetries are in general broken in the superfluid states, and symmetry 
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considerations [l8|, (and references therein) have been used to find the inert states. 
These considerations have also been applied to "unconventional" superconductivity in solid 
state systems, where the superconducting state has symmetries lower than the underlying 



20. 



2l| . In this case, the normal 



crystals, in particular in the context of heavy fermions, e.g. 
state possesses only discrete rather than continuous rotational symmetries. 

In this paper, we shall carry out such a program for the spin condensates, finding the inert 
states and their symmetries. We shall limit ourselves to the case of no external magnetic 
field (linear and/or quadratic Zeeman fields). We discuss the general theory behind this 
work in Sec [Til Results up to S = 4 are then presented in Sec IIII1 We then apply these 
results to discuss vortices on Sec [TV] Sec [V] provides a brief summary, and some details are 
relegated to the appendices. 

II. BASIC THEORY 

We here outline the basic theory behind this work. The present procedure of finding the 



inert states is based on a theorem by Michel 22J. This theorem has been discussed and 



explained in, e.g., 19|. Rather than restating this theorem in its full generality in abstract 
terms, we shall simply describe it physically for the case of S = 1 here. Since S = 1 has 
the same symmetry as angular momentum L — 1, we shall use these two terminologies 
interchangably. In the uncondensed phase (the normal phase), the system is invariant under 
all spatial rotation SO(3) and gauge transformation U(l), and time reversal O. (There is no 
need to consider inversion since we shall always consider fixed L, so all states have definite 
inversion symmetry). We are interested in finding the stationary points of a free energy 
F[^f] which is a functional of F[^f] is invariant under any transformation ^ — > g^> where 
g is any operation in the group G = SO (3) x U(l) x ©. For L = 1, it can be shown easily 
that all order parameters can be written in the form k ■ (rh + ch)e n where rh and n are real 
orthogonal unit vectors, and c is a complex number with |c| < 1. (We shall, for simplicity of 
presentation, not always normalized our wavefunction, since this normalization has no effect 
on symmetries. We can restrict ourselves to |c| < 1 because if otherwise, we can always 
divide the state by c and then rewrite the state in the above form). It is also sufficient to 
consider only the two cases c = and Imc ^ (if Imc = 0, then the state is just (k ■ m')e n 
for some rh'). Special choices of rh, n and phase factor 7 results in reference states, which 
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we shall take to be k ■ z and k • (x + cy). All states (k ■ m)e' n can be obtained from the 
reference state k ■ z by rotation and gauge transformation. Equivalent statement applies to 
the states k ■ (m + cn)e n for a given c. [Technically, the collection of states {g^} obtained 
from \1/ by operating the transformation g on it are said to form the orbit of ty] . As far as our 
problem of obtaining stationary states are concerned, these solutions related by symmetries 
are "equivalent" solutions. ^ is an inert state if and only if g^/ is an inert state. 

An important consideration in determining whether a particular state \l/ is an inert state 
is the symmetry group under which ^ is invariant, that is, the isotropy subgroup H of 
where H contains all operations h in G such that h^> = ^ . We first note here that a state 
which can be obtained by a transformation g on f has symmetry simply related to that of 
the original reference state ^> . [Technically, g^ has the isotropy group gHg~ l thus these 
isotropy groups are said to be conjugate to each other]. Consider now the reference state 

= k • (x + cy). If c = ci + ic2 with c\ 7^ and c 2 7^ then H contains the identity E alone. 
All states with c\ 7^ and c 2 7^ have the same symmetries. [Technically, the collection of 
orbits {g^ c } from states \l/ c with the same symmetry are said to form a stratum}. For any 
given (01,02) 7^ (0,0), we can choose a (c'^c^) arbitrary close to the original (01,02) such 
that (7^ g^ c for any g in G) has the same symmetry as the original ^> c . [Technically, the 
orbit {g^ c } therefore not isolated in its stratum] . Michel's theorem then says that the state 
\1> C is not an inert state. For example, we do not expect the state k • (x + ^y) with c = ^ 
to be an inert state. While it may minimize a particular free energy F, there is no reason 
to expect c to remain at the same value of c = ^4= when the parameters in the free energy 
F are changed. Similar remarks apply to the case where c = ic2 with | C2 j < 1 (though H 
contains also other elements in additional to E). 

The states with, however, c = i is special in that the isotropy group H would then become 
{R z (a)e~ ia , £/ 2 7 e~ 2i7 6} where R z (a) denotes a rotation of angle a about z, and U% denotes a 
7r rotation about the horizontal axis 7 = cos7x + sin7y. (Note that 9 does not commute with 
e~ l27 . We are using the notation where operations start from the left). A small modification 
of c away from % would result in a state that has a completely different symmetry (and 
is not g^ for some g in G, that is, not simply another order parameter related to \l/ by 
a transformation in G). [Technically, we say that the orbit g^/ is isolated in its stratum}. 
Michel's theorem then states that ^> c =i is an inert state. This state is a stationary point of 
the free energy irrespective of the parameters that enter in F. 
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Similar statement applies to \I/ = k ■ z, (related to the state ^ c =o by a rotation) where 
the isotropy group is H = {R z (a), U^e 17 *} x O. (From the above argument, we have in fact 
shown that the only inert states for L=l are the axial state k ■ (x + iy) and the polar state 
k ■ z. In the spin wavefunction language they are (1, 0, 0) and (0, 1, 0) respectively with 
isotropy groups already given above). 

Therefore, in general (not restricted to L = 1), one is searching for states \I/ that have 
special symmetries in the sense that a slight deviation of the coefficients entering \l/ would 
give another state with an isotropy group different from H (but barring cases where 
is just an "equivalent" solution, that is, ^' = and thus just a symmetry-transformed 

i.e., cases where the new order parameter is in the same orbit and thus with symmetry 
group gHg -1 ) [Technically, we look for states \l/ such that its orbit is isolated in its stratum}. 

The procedure is thus to search for all subgroups H of the combined group SO(3) x U(l) 
xG with broken U(l) gauge symmetry which would determine \I/ uniquely. Each of these \I> 
would then be an inert state. (In principle it is sufficient that H determines a set of possible 
\l/'s, such that the orbits of these order parameters are isolated. However, in the present 
case of interest, this ambiguity is only in an overall ± sign for the wavefunction. Thus 
effectively each of these H determines a unique inert state \1/.) It is useful to note that once 
such an H is found, then there is no need to consider the subgroups of H. This is because 
any ^2 invariant under a subgroup H2 of Hi but not Hi itself cannot yield an inert state. 
(To see this, consider \I/(A) = ^2 + A\&i where ^1 (7^ ^2) is invariant under Hi. Then ^(A) 
is invariant under H2 for any A and hence, by the argument given, ^2 cannot be an inert 
state.) 

The subgroups of SO(3]x U(l) x6 that break the gauge symmetry can be constructed 
in a similar manner as in |20|]. There, Volovik and Gorkov constructed all possible isotropy 
subgroups for unconventional superconductors where the normal state possesses discrete 
rotational symmetries of crystals instead of SO (3). A principal observation, which still 
applies here, is that U(l) is isomorphic to rotation about a fixed axis. Thus the required 
groups would involve elements that are in general product of a spatial rotation and a gauge 
transformation. 

In our present case, these subgroups fall into two categories. The first involves continuous 
groups. These can be constructed easily (see below). The second involves discrete groups. 
Our task is simplified since the discrete subgroups of SO (3) are well known. They are the 
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three groups of regular polyhedrons: O of the octahedron, T the tetrahedron, and Y the 
icosahedron; the dihedral groups D m and the cyclic groups C m . 

Before we actually list these groups we note further a great simplication. Let us consider 
a general spin wavefunction \I/ = (£5, C-s)- If the group H contains the symmetry element 
h = C n e tj2n / n where j is an integer which can be restricted to — n < j < n (j is an integer 
follows from the requirement that h n must also be in H: as C™ = 1, we must have j being an 
integer), since under C n , ( m — ► Cme~ im2n ^ n , we must have e l U- m ) 27r / n = 1 f or £ m ^ 0. Hence 
only (j, Cj±n e ^c can be finite. If there is only one such finite ( m then the order parameter 
has been determined uniquely. The state then has in fact continuous symmetry R z (a)e^ a . 
If there are more than one non-zero ( m , say Q and Q- n , then it can be an inert state only if 
there are further other symmetry elements in H that puts a constraint relating Q and Q-n- 
Only when j — n = j can this be achieved by a horizontal rotational axis, and the symmetry 
group would then be a combined group constructed from the dihedral groups D n (or a larger 
group containing this). If j — n ^ j, such requirements would arise only from groups that 
are constructed from the regular-polyhedral groups O, T and Y. 

Hence we can follow the following procedure. First we consider all the continuous isotropy 
subgroups and find the corresponding order parameters. Then we consider the (non-dihedral) 
subgroups that are constructed from the polyhedral groups: O x U(l) x6, T x U(l) X0, 
and Y x U(l) X0. If an order parameter for a particular isotropy group can be found, then 
we discard all its subgroups. If not, then we go to a group of lower symmetry. This is basically 
the only part of the procedure that is more involved. After such groups are exhausted, the 
only inert states that can be left are those correspond to the dihedral symmetries. These 
states can be obtained by inspection (with the procedure already hinted above in the last 
paragraph: see examples below). For our case with L < 4, we shall argue in Appendix 
[A] that five-fold symmetry elements cannot be involved, hence there is no need to consider 
the icosahedral group Y (and also dihedral groups containing C5). The only non-dihedral 
polyhedral groups needed are thus those constructed from O x U(l) xO and T x U(l) x0, 
and they can already be found in [20]. 

Now we proceed therefore to list the continuous subgroups, and subgroups which can be 
constructed from the polyhedral groups O and T. The continuous ones are easy. We can 
combine R z (a), continuous rotation of angle a about the axis z, with gauge transformation 
e %na w j iere n j s an integer. The only state that would be invariant under R z (a)e ma for all 
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a is the wavefunction where the only non-zero component ( m in \1/ is ( n . Corresponding 
wavefunction of angular momentum L is ( n Y£(k). This wavefunction automatically has 
symmetry elements around horizontal axes 7 combined with suitable phase factors (due to 
the definition from spherical harmonics). For L — n even, the isotropy group is 

H* = {R z (a)e ma , £/ 2 V iri7 0} (1) 

whereas for L — n odd, we have 

H° = {R z (a)e ma , e l(2n ^ +n) Q} (2) 

We need only consider non-negative ra's, since the state with only (_ n being non-zero 
is related to the state with ( n being the only non-zero element via either the time-reversal 
symmetry or a horizontal rotation by 7r. 

For the discrete (non-dihedral) subgroups contructed from O and T, they have already 



been considered by Volovik and Gorkov |20| and so we simply list their answers below: 



0x0 = {E, 8C 3 , 3C 2 , 6U 2 , 6C 4 } x (3) 

This group corresponds to just breaking the gauge symmetry in O x U(l) xO. Here C3 are 
three-fold rotations along the diagonals, U2 are two-fold rotations along the edges, and C4 
are four-fold rotations about the faces of a cube, and C2 = C\. 

0(T) x = {E, 8C 3 , 3C 2 , 6U 2 e l7r , 6C 4 e l7r } x (4) 
This group is isomorphic to O, but contains only T as a subgroup. 



0(D 2 ) = {E, 3C 2 , 2Q'0, 2C*|e0, 2Qe 2 0, 

4C 3 e, 4C|e 2 , 2U^Q, 2U^ y e&, 2l^e 2 0} (5) 

This group contains D 2 as a subgroup. Here U 2 denotes a ir rotation about the edge of 
a cube perpendicular to x (i.e., x — > —x and y <-> z or y <-» —z), and e = e 2t7T ^ 3 . 
We also have the tetrahedral groups 

T x = {E, 3C 2 , 4C 3 , 4C 3 2 } x (6) 
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and 



T(D 2 ) = {E,3C 2 ,4C 3 e,4Cle 2 } (7) 

This last group is isomorphic to T but only contains D2 as a subgroup. We note here also 
that T is a subgroup of both O and O(T) and T(D 2 ) is a subgroup of 0(D 2 ). It turns out 
that neither T nor T(D 2 ) give rise to inert states for any of the cases we have examined 
below. 

This complete the list of all non-dihedral groups that we would need (for L < 4). 



III. RESULTS 



Now we apply our theory to, in order S (or L) = 2,3,4. For simplicity in notations, 
when writing the wavefunction for L, we would simply write x, y, z instead of k x , k y , k z . It 
is understood that x 2 + y 2 + z 2 = 1. 



A. 5 = 2 



The inert states for L = 2 corresponding to d-wave Cooper pairing has in fact been 



investigated before by [18 



(see also 



191 ]) by considering the system as strongly spin-orbit 



coupled L = 1 and 5=1 Cooper pairs with total angular momentum J = 2. Here we would 
like to present our arguments directly in terms of the procedure outline above so that it can 
be easily compared with the results for S = 3 and 4. 

Continuous Groups: We can easily write down the three order parameters corresponding 



to continuous groups H^'° for n = 2, 1, 0. We shall call these states 
indicates ferromagnetic and P polar, using the language as in, e.g., 
state wavefunctions are 



F 1 and P (Here F 



23] ). Explicitly, the 



F 2 = (1,0,0,0,0) 
F x = (0,1,0,0,0) 

P = (0,0,1,0,0) (8) 
with corresponding polynomial forms (x + iy) 2 , —z{x + iy) and (2z 2 — x 2 — y 2 ). 

8 



Discrete groups: It can be easily verified that no order parameter can satisfy the groups 

Ox@, 0(T) x9 and T x6. To see this, we note that all L = 2 wavefunctions can be 

written as polynomials of second degree in x, y, z. The elements 3C 2 require \l/ to be a linear 

combination of x 2 , y 2 , z 2 only, (e.g., xy is odd under C 2 and hence must be rejected). All 

these three groups contain eight three-fold rotations C3 along the diagonal of a cube, where 

(x, y,z) — > (y,z,x), (x,y,z) — > (—z, —x,y) etc. The only function that has this property is 

x 2 + y 2 + z 2 , which however does not belong to L — 2. 

The discrete group 0(D 2 ) in eq (El) uniquely determines (up to a ± sign) the order 

I 

parameter to be x 2 + ey 2 + e 2 z 2 (see also [20|). We have just seen that 3C 2 limit us to the 
polynomials x 2 , y 2 and z 2 . Under C3, we have the mapping (xyz) — > (yzx) and hence the 
elements C^e requires that the wavefunction be proportional to x 2 + ey 2 + e 2 z 2 . Elements 
such as C46 fixes the overall phase factor to be real (under Cf, (x,y,z) — > (x,z, —y))- We 



can also check that elements such as U~2® are satisfied. As before 



23|, we shall refer 



to this state as the "cyclic" state. We shall denote it by C, and further discuss this state 
below. 

Since we have already found the order parameter for 0(D2), there is no need to consider 
its subgroup T(D 2 ). This completes the non-dihedral polyhedral groups. 

It remains then to consider the ones corresponding to order parameters where only (± m 7^ 
for some m^O and symmetries related to the dihedral groups. For m = 2 we can have 
the state 

P 2 = (1,0, 0,0,1) (9) 
Its polynomial form is x 2 — y 2 and easily seen to possess the symmetry (same notation as 



111 



20|) 



D 4 (D 2 ) x 6 = {E, C 2 , 2U 2 , C 4 e m , 2U' 2 e™} x 6 (10) 

It can easily be checked that this group actually determines the order parameter uniquely 
(again up to ± sign, a warning that we would not repeat again below) to be x 2 — y 2 , and 
hence the state is inert. Note that a rotational about z changes the relative phase between 
(±2, thus all states of the form e l ^(e~ 2ta , 0, 0, 0, e 2ta ) are related to P 2 by operations within 
G. States with |£_ 2 | 7^ IC2I as the only non-zero elements have lower symmetries (in the sense 
of subgroup) and therefore need not be considered. 
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One may attempt to try the state \I/ = (0, 1, 0, 1, 0), which in polynomial form reduces to 
—i2zy. Thus up to a gauge tranformation and a rotation, this state is just the same as 2xy. 
We note here that a 7r/4 rotation about z changes x 2 — y 2 to 2xy. Thus the state (0, 1, 0, 1, 0) 
can be transformed to the state P2 and does not correspond to a new inert state. 

This completes our search for inert states for L = 2. We found five inert states: F 2 , F\, 
Po with continuous symmetries, C and P2 with discrete symmetries. 

In previous wo^aCi,^— ion of^a^e free energy were per- 
formed. Our inert states accounted for all minina found in 0, [2]. There, it turns out 
that Fi is never the absolute minimum for any parameters entering F"[\I/]. It also turns out 
to be a peculiar fact that, to this order, P and P2 are always degenerate. Furthermore, 



any real linear combination between P and P2 are also degenerate [15| [Mermin has shown 
that, when the state is written in the form B^XiXj where £1,2,3 = x, y, z, any real symmetric 
traceless B^ produces the same free energy. Since any such matrix B^ can be diagonalized, 
we see that these \I/ can always be recast into the form a(2z 2 — x 2 — y 2 ) + b(x 2 — y 2 ) with 
a and b real. This is thus a real linear combination of Po and P2] We however expect that 
these accidental degeneracies would be lifted in more general free energies (e.g. containing 
higher order terms in 

Let us further comment on the cyclic state. It can easily be checked that the polynomial 
ex 2 + e 2 y 2 + z 2 is proportional to V2Y® + i(Y 2 + Y^ 2 ), and hence the state vector can be 
written as 

C = (i, 0, V2, 0, i) (11) 

corresponding to a form used in [8j. A 7r/4 rotation around z can transform this into 
(1,0, y/2,0, -1), that is y/2Y£ +JY 2 - Kr 2 ). We remark here that this state was in fact 



found by Anderson and Morel 



24j to be the minimum energy state for a d-wave Fermi 



superfluid within weak-coupling BCS theory. A pictorial way of showing the symmetry of 



this state can also be found there (see also 



20|). 



Since the group 0(D 2 ) has explicitly the symmetry element C^e, we know immediately 
that, if one uses quantization axis along this axis, the wavefunction must assume the form 
(0, Ci, 0, 0, C-2), z.e., only d and C_ 2 are non-zero. Indeed, let us take w = (x + y + z) / \f?> and 
u = (— x + y)/y/2, and v = w x u to be the new orthonormal axes. Using the transformation 
from r = xx + yy + zz to the new coordinates f — uu + vv + ww, simple algebra shows 
that ex 2 + e 2 y 2 + z 2 = —\[2%w{u + iv) — — ™^ ■ Therefore, with w as the quantization 
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axis, the wavefunction can then be also written as (0, v2i, 0, 0, —1), the form as advertized. 
Of course, there are also many other possible forms for C (see., e.g., [h]]). We remind the 
reader here that this state is unique up to gauge and rotations (as also proven by Mermin 

Q)- 

B. 5 = 3 

The states that have continuous symmetries are 

F 3 = (1,0,0,0,0,0,0) 
F 2 = (0,1,0,0,0,0,0) 

Fi ee (0,0,1,0,0,0,0) (12) 
P ee (0,0,0,1,0,0,0) 

with corresponding polynomial forms — (x+iy) 3 , z(x+iy) 2 , — (5z 2 — l)(x+iy), and (5z 3 — 3z). 

To construct the states belonging to the polyhedral groups, it is convenient to first notice 
that an L = 3 state can always be written in the form of a cubic polynomial in x, y, z, 
thus Bij k XiXjXk where £1,2,3 = x, y, z. We can further limit ourselves to B^ k symmetric 
under all possible interchanges of indices ijk. Orthogonality to L — 1 then requires < 
BijkXiXjXkXi >= where < ... > indicates an angular average. Hence we have Ylii^m = 
and so it is sufficient to specify only seven (=2x3+1) coefficients, which can be chosen 
to be B xyz , B xyy , B xzz , B yxx , B yzz , B zxx , B zyy (by eliminating e.g., B xxx in terms of B xyy 
and B xzz .) All the polyhedral groups OxO, O(T) xO, 0(D 2 ), T x and T(D 2 ), contain 
3C 2 . The only polynomial that satisfies this is xyz. (e.g., invariance under C z leaves only 
xyz, xyy and xzz. The latter two are eliminated by the requirement of invariance under 
Cf.) Therefore obviously O, 0(D 2 ), T(D 2 ) cannot be satisfied. However, it does satisfy 
O(T) xO. We can check in fact O(T) xO uniquely determines the order parameter to be 
xyz. Note that xyz = zlm(x + iy) 2 , hence the spin-wavefunction for xyz can be written as 
(0, — i, 0, 0, 0, i, 0). Though by a suitable rotation about z, we can change this to the form 
(0,1,0,0,0,1,0) and thus call this P 2 , for reasons to be explained below we shall instead 
call this state D: 
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D = (0, —i, 0, 0, 0, i, 0) (13) 

Since T x9 is a subgroup of O(T) x6, we have exhausted all states under the non- 
dihedral polyhedral groups. Hence we proceed to the states with only (± m being finite 
obeying dihedral groups. For m = 3 we can have the state (—1, 0, 0, 0, 0, 0, 1), and we shall 
call this P3. It has the polynomial form Re(x + iy) 3 and the isotropy group 

D 6 (D 3 ) x = {E, 2C 3l 3U 2 , C 2 e in , 2C 6 e*«, x 6 . (14) 

We can check again that D 6 (D 3 ) uniquely determines the order parameter to be Re(x + iy) 3 . 
For m = 2, we would obtain the state P2 which was already discussed under the group O(T) 
xG. For m=l, the wavefunction (0, 0, 0, 1, 0, 1, 0, 0, 0) would then have the form (5z 2 — l)y 
up to a phase factor i. This state is equivalent to (5z 2 — l)x which has the symmetry 
D2(U 2 ) x 6 = {E, C 2 e l7T , U 2 , U 2 e l7T } x 6 and is thus a subgroup of D 6 (D 3 ) x 6 considered 
before and thus this state is not an inert state. We thus have exhausted the list for states 
with dihedral symmetries. 

Summarizing, for S = 3 we found six inert states, F 3 , F 2 , F x and Po with continuous 
symmetries, D and P 3 with discrete symmetries. 



Minimization of a general quartic free energy has been carried out by [12j, ll3|] • m > 
F 3 , F2 and P 3 appeared explicitly as absolute minima for suitable interaction parameters. 
They also find a state (which they named state D) which we shall call D' with the spin 
wavefunction (y/2, 0, 0, iy5, 0, 0, y/2). The form suggests that it is an inert state. Since it 
has three-fold symmetry but it is obviously not P 3 , this suggests that it should be the state 
D expressed in rotated axis. We shall verify that this is indeed the case in Appendix Q31 
Thus their state D' has isotropy group (conjugate to) O(T) X0. 
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C. 5 = 4 



For continuous groups, we have 

F 4 = (1,0,0,0,0,0,0,0,0) 
F 3 = (0,1,0,0,0,0,0,0,0) 

F 2 = (0,0,1,0,0,0,0,0,0) (15) 
Fx = (0,0,0,1,0,0,0,0,0) 
P = (0,0,0,0,1,0,0,0,0) 

(16) 

with polynomial forms (x + iy) 4 , —z(x + iy) 3 , (7z 2 — l)(x + iy) 2 , — (7z 3 — 3z)(x + iy), 
(35,2 4 — 30z 2 + 3). For discrete groups, we have to consider even order polynomials up to 
fourth order (we shall require orthogonality to L = and 2 later). For the polyhedral groups 

xO, O(T) x0, 0(D 2 ), T xG and T(D 2 ), the elements 3C 2 eliminated all polynomials 
except x 4 , y 4 , z 4 , y 2 z 2 , z 2 x 2 , x 2 y 2 , x 2 , y 2 , z 2 . To obey 8C3 in O x9 and O(T) x0, we must 
then have the linear combination 

= c (x 4 + y 4 + z 4 ) + a(y 2 z 2 + z 2 x 2 + x 2 y 2 ) + b(x 2 + y 2 + z 2 ) (17) 

Obviously O(T) x6 is not satisfied (no required sign changes under C 4 ), but O x 6 is if c, a, b 
are all real. Before we continue, we note that the three functions multiplying the coefficients 
c, a, b are not linearly independent. Indeed, using x 2 + y 2 + z 2 = 1 we can verify that 
(x 4 + y 4 + z 4 ) = (x 2 + y 2 + z 2 )—2(y 2 z 2 + z 2 x 2 + x 2 y 2 ). Let us choose to use (y 2 z 2 + z 2 x 2 + x 2 y 2 ) 
and x 2 + y 2 + z 2 = 1 to be the two independent functions, thus setting c = in eq ffTTj) . 
Orthogonality to L = and 2 requires < ^ >= , < ^x 2 >=< ^y 2 >=< ^z 2 >= 0. 
Evaluating the angular averages 29(, we find that all these relations are satisfied if b = —a/5. 
(The fact that these equations are dependent is not a coincidence, since < x 2 {x 4 + y A + z A ) > 

1 < (x 4 + y 4 + z 4 ) >= < x 2 (y 2 z 2 + z 2 x 2 + x 2 y 2 ) > / < (y 2 z 2 + z 2 x 2 + x 2 y 2 ) > =< x 2 >= 1/3 
etc from symmetry.) Thus OxO determines the order parameter uniquely and the state is 
inert. Choosing a = 5 gives us the wavef unction 

= 5( y V + z 2 x 2 + x V) - 1 (18) 
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For lack of a better name we shall just call this state 0. In Appendix Owe shall verify that 
the corresponding spin wavefunction, up to a normalization factor and sign, is given by 

O = (V5, 0, 0, 0, Vli, 0, 0, 0, \/5) . (19) 

We need not consider T x since it is a subgroup of O x 0. 

In 0(D2) or T(D2), the element C%e require that we have instead the linear combination 

$ = c(x 4 + ey 4 + eV) 

+a{y 2 z 2 + ez 2 x 2 + e 2 x 2 y 2 ) + b{x 2 + ey 2 + eV) (20) 

So far, c, a, b can be complex. This state obviously satisfies T(D 2 ) (after orthogonality to 
L = and 2 have been shown). To satisfy the larger group 0(D 2 ), C%Q requires that 
the coefficients c, a, b are all real. [In this case, one can check that elements such as U^ X Q 
(x <-> —x, y <-> z and then complex conjugate) are also satisfied.] Before we proceed, we 
note again that the three functions multiplying c, a, b in eq ( 1201) are not linearly independent. 
Indeed, using x 2 + y 2 + z 2 = 1 and 1 + e + e 2 = 0, we can verify that (x 4 + ey 4 + e 2 z 4 ) = 
(x 2 + ey 2 + e 2 z 2 ) + (y 2 z 2 + ez 2 x 2 + e 2 x 2 y 2 ). Let us proceed to work with the form with c = 0. 
It remains to find a, b to satisfy orthogonality to L = and 2. < ^ >= is trivial (since 
1 + e + e 2 = 0). < ^x 2 >= and < % 2 >= both require -a + 7b = (< ^z 2 >= 
is then automatically satisfied), and therefore both a and b can be chosen real. Hence we 
obtain 

* = 7(y 2 z 2 + ez 2 x 2 + e 2 x 2 y 2 ) + {x 2 + ey 2 + eV) (21) 

as a state satisfying 0(D 2 ) (thus there is no need to consider T(D 2 )). In appendix O, 
we show that, up to a renormalization constant and phase factor, the corresponding spin 
wavefunction is 

C = (V7, 0, y/Ui, 0, -VW, 0, Vl2i, 0, y/7) (22) 

This is our state with 0(D 2 ) symmetry. We are calling this state cyclic (C) in analogy with 
the cyclic state in L = 2. This completes our consideration of polyhedral groups. 

The states with dihedral groups can be found as in L = 2 and 3. The state P4 = 
(1, 0, 0, 0, 0, 0, 0, 0, 1) has polynomial form Re(x + iy) 4 and symmetry 

D 8 (D 4 ) x = {E, C 2 , 2C 4 , 4C 8 e^, 4£/ 2 , AU'^} x (23) 
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The state P3 = (0, i, 0, 0, 0, 0, 0, —i, 0), (related to (0, 1, 0, 0, 0, 0, 0, 1, 0) by a rotation around 
z), polynomial form zlm(x + iy) 3 , has symmetry D@(D3) x (eq (j!4p ). We need not consider 
(0, 0, 1, 0, 0, 0, 1, 0, 0) since then it has the form (7z 2 — l)(x 2 — y 2 ) and symmetry D 4 (D 2 ) x 6, 
and thus is only a subgroup of D 8 (D 4 ) x O. Also, for the similar reason there is no need 
to consider (0, 0, 0, 1, 0, 1, 0, 0, 0) which is proportional to (7z 3 — 3z)y which has symmetry 
D 2 (C2) x 6. This completes the list of discrete groups. 

Summarizing, we found nine inert states: F4, F3, F2, Fi, Po with continuous symmetries 
and O, C, P 4 , P 3 with discrete symmetries. 



IV. VORTICES 



As an application of our results, we consider topological line defects (vortices) of the spin 
condensates. The distinct types of vortices have already been considered in |6| for S = 1 
and for the states F 2 , Fi and C of S = 2. 

The possible types of vortices are determined by the map of a circle in real space to 
the space of non-equivalent order parameters. We apply the theorem that has been proven 



in, e.g, [271 ]. which states that 7ri(G'/H') = H'/H' . Here, G' is a simply-connected group 
characterizing the symmetry of the fully symmetric (normal) phase, and H' is the isotropy 
subgroup (within G') of the state under consideration, and H' is the part of H 7 that is 
connected (in the sense of a neighborhood) with the identity E. If H'/Hq is abelian, the 
line defects are simply classified by the elements of this group. Otherwise, the topologically 
distinct type of vortices are classified by the conjugate classes of this group. 



A. S=2 



This has already been treated in 



Q (c.f. 



17j) except the states P 2 and P . For ease of 



comparison with S = 3 and 4 however we shall also repeat the results for the other states 
below. 

For the states Fi 2 , since a rotation about the z axis would also generate a gauge trans- 
formation, it is simpliest 27j to just consider the symmetry group SO (3) for the normal 
state. For Fi, a reference state can be taken just as z(x + iy). Then H contains the identity 
E alone. To have a simply-connected G', we "lift" G = SO(3) to G' = SU(2), and then H' 
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becomes H' = {E, Q} where Q ^ E and Q 2 = E. H consists of E alone and H'/H = H'. 
We thus have two classes of vortices for Fi. 

For F 2 , we take the reference state (x + iy) 2 . Then H = {E,R z (tc)} consisting of two 
classes. When we lift from G = SO (3) to G' = SU(2), each element in H is doubled. Thus 
there are four classes of vortices. 

For Pq, we choose the reference state (2z 2 — x 2 — y 2 ) and take G = SO(3) x U(l). Then 
H = Hq = {R z (a),U 2 } (see eq (JTJ consisting of two classes. We need however a simply- 
connected G', which can be taken to be G' = SU(2) x T^, where is the translational 
group for the phase 0. Then H' = {u(z, a)T 2nn , u(j, ±n)T 2nn }. Here, we have used the 
representation u(n, a) = cosf +in-<rsm% in terms of Pauli matrices for the rotations in SU(2). 
Hq = {u(z, a)}, and thus H'/H consists of the elements (cosets) {u(z, a)T_ 27r }, {u(z, a)}, 
{u(z, a)T +2n }, {u(j, ±7r)T 27r }, .... Since a product between u(z,ai) and u(z,a 2 ) or 
between u(7i,±7r) and m(7 2 ,±7t) both yield u(z, a) for some a; and the product between 
u{z) and u{fy, ±7r) yield another element ^(7', T 71 "); H'/H is isomorphic to Z 2 x Z, where we 
can assign the integer to {u(z, a)} and 1 to {tt(x for the group Z 2 and the circulation 
number n for the group Z. Note this is different 16( from the polar state Po for spin 5=1. 
The source of difference is that the isotropy group for the state P is Hq for S = 2 but Hq 
for 5 = 1. 

The cyclic state C has been treated in [if]] using the language of the spin- wavef unction 
and its rotational properties. Here the symmetries we obtained before becomes handy. 
If G = SO(3) x U(l) x 6, we have already found H = 0(D 2 ) (eq To consider the 

vortices, we only need to consider instead G = SO(3) x U(l), and then we obtain H = 
{E, 3C 2 , 4C 3 e, 4C|e 2 }. This group is isomorphic to the tetrahedral group with T = 
{E, 3C 2 , 4C 3 , 4C|} which has four conjugate classes. The double group T' corresponding 



T has seven classes 



28|, {E},{Q},{3C 2 ,3C 2 Q},{4C 3 },{AC 3 Q},{4C 2 },{4C 2 Q}, as each 



class in T is doubled except the one ({3C 2 }) which consists of bilateral two-fold axes [26]. 
When G = SO(3) x U(l) is lifted to the G' = SU(2) x T , we get H' = {E}, H'/H = H' 
where the group H' has classes {ET 2n7T }, {QT 2nn }, {3C 2 T 2nn ,3C 2 QT 2n7T }, {4C 3 T (2n+ 2 )7r }, 
{AC 3 QT [2n+ 2 )7T }, {4C|T (2n+ 4 )7r }, {4CfQT (2n+ 4 )7r }, where n is an integer. Thus the vortices 
are divided into seven types in additional to the circulation number n. Note however due to 
the presence of e = e t2n ^ 3 factors in H and hence H', H'/H is not isomorphic to the product 
of the double group T' times the discrete translational group T 2niT . 
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The polar state P2 with discrete symmetry can be treated similarly. If 
G = SO(3) x U(l), then H is simply D 4 (D 2 ) = {E, C 2 , 2U 2 , 2C 4 e i7r , 2U' 2 ^} 
isomorphic to D 4 = {E, 2C 2 , 2U 2 , 2C 4 , 2U 2 } with five classes. The 
double group D' 4 corresponding to D 4 consists of seven classes 28]: 
{E}, {Q}, {C 2 , C 2 Q}, {2U 2 , 2U 2 Q}, {C 4 , C 3 Q}, {Cf , C 4 Q}, {2U' 2 , 2U' 2 Q}. When we lift 
to the connected group G' = SU(2) x T^, we obtain H = {E}, and H'/H = H' consists 
of the following classes: {ET 2nn }, {QT 2nn }, {C 2 T 2n7T , C 2 QT 2n7T }, {2U 2 T 2n7T , 2U 2 QT 2n7T } , 
{C 4 T 2n7r , ClQT 2n7T } } {C|T(2 n +i)7r, C^QT^n+i)-*} > {2U' 2 T( 2n+ i^ 1T ,2U' 2 QT{ 2ri+ i^ 1T }. There are 
seven vortex types in additional to the circulation numbers n. 

B. 5 = 3 

For the state F 3 , it is again simpliest just to use SO(3) and the reference state (x + iy) 3 . 
H consists of three elements, E,R z (2Tr/3),R z (4ir/3). The corresponding double group H' 
has six elements. There are six distinct classes of vortices. 

F 2 and Fi has the same corresponding H as their S = 2 counterpart. The distinct classes 
of vortices are again four and two respectively. 

For the polar state Po, we choose the reference state (5z 3 — 3z) and take G = SO(3) 
x U(l). Then H = {R z (a),U 2 e ln } (see eq 02)) consisting of two classes. When using 
the connected G' = SU(2) x T^, H' = {u(z, a)T 2n7r , u(j, ±7r)T( 2n +i)7r}- H' = {u(z, a)}, and 
thus H'/Hq consists of the elements (cosets) {u(z, a)T_ 27r }, {u(z, a)}, {u(z, a)T +2n }, 
{m(7, ±7r)T^}, Note now that the cosets with 7/(7, ±7r) are associated with translation of 
phase of odd integral multiple of ir. Again noting the properties of the products of u(n, a) 
discussed in S — 2, we find that H'/Hq is isomorphic to Z, as it is the case [3] for the state 
P belong to S = 1. 

For the state D, when we just use G = SO(3) x U(l) then H = O(T) = 
{E,8C 3 ,3C 2 ,6U 2 e in ,6C 4 e ln } is isomorphic to O = {E, 8C3, 3C 2 , 6U 2 , 6C 4 } with thus 
five classes. The double group 0' corresponding to O has eight classes, (only {E} , 
{8C3} and {6C 4 } are doubled). Vortices can therefore be classified into the follow- 
ing classes: {ET 2nn }, {QT 2n7T }, {4C 3 , 4CfQ}T 2n7r , {4C|, 4C 3 Q}T 2n7r , {3C 2 ,3C 2 Q}T 2n7I , 
{3C 4 , 3C|Q}T (2n+1)7r , {3C|, 3C 4 Q}T ( 2n+i)-x, {6U 2 ,6U 2 Q}T {2n+1)n . There are thus eight 
classes in additional to the circulation number n, and those elements originally associated 
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with e m in the group O(T) are now associated with phase translations (2n + 1)tt. 

For the polar state P3, if we just use G = SO(3) x U(l) then H = D 6 (D 3 ) = 
{E, C 2 , 2C 3 , 2C 6 e" r , 3U 2 , 3U 2 e ln } was already given in eq (III]) . This group is isomorphic 
to D 6 = {E, C 2 , 2C 3 , 2C 6 , 3U 2 , 3U' 2 } with six classes. The double group D' 6 corresponding 
to D 6 has nine classes: all classes are doubled except C 2 , 3U 2 and 3^2 • Vortices are there- 
fore divided to nine classes, in additional to an integer n specifying the circulation. (Those 
classes that correspond to C§e m and U^e™ are associated with phase translations T( 2n +i) 7r 
rather than T 2n7r (c.f. the state D). 

C. S = 4 

These can be obtained similarly and so we shall be brief. Fi, F2, F3, F4 have respectively 
two, four, six and eight classes of vortices. Po has symmetry same as Po of S = 2 and thus 
with the same classes of vortices. A parallel statement applies to the cyclic state C. The state 
O has eight classes of vortices in additional to the circulation numbers. We note however 
that all vortices for O have winding phase factors 2nir, thus different from the state D of 
S — 3 where the symmetry group is O(T) instead of O. The state P4 has symmetry group 
D 8 (D 4 ). The double group D' 8 has nine classes. So vortices for P 4 has nine classes of vortices 
in additional to the circulation number, and each class is associated phase translations of 
(2n + 1)tt or 2mr depending on whether the factor e l7T appears or not in Dg(D4). The state 
P 3 is analogus to that in S = 3. 

V. CONCLUSIONS 

To conclude, we have listed all inert states for a spin superffuid for S < 4. These states 
are always stationary points of any general free energies. We have also given the residual 
symmetries (hence broken symmetries) of these inert states. 

This research is supported by the National Science Council of Taiwan under Grant No. 
NSC95-2112-M-001-054. 
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APPENDIX A: GROUPS WITH FIVE-FOLD SYMMETRY ELEMENTS 



Here we explain why groups containing five-fold symmetry elements C5 (without contain- 
ing the continuous symmetry Coo), i n particular the icosahedral group Y, are not involved 
for inert states with L < 4. First, it is clear that, since L < 4, the wavefunction cannot be 
invariant under C5 itself. Next we show that H cannot contain G itself. For if it did, then it 
would contain both C^e*^ 5 and C^e^Q for some The latter element can be rewritten as 
OCse - ^ 5 and hence H would contain Cf , leading to a contradiction. Similarly we can show 
that H cannot contain an element of the form h = C^e^Q. If h is an element in H, then 
so must h 5 = e^ 5 G. Thus apart from a gauge transformation, the state \I/ is time-reversal 
invariant = v|/e^ 5//2 is invariant under 0). which is again not possible. Therefore C5 
must be in combination with a gauge transformation only. As seen before, this phase factor 
must be integral multiple of 2n/5. Hence let us consider the element C^e*- 7271 "/ 5 where j is 
an integer. As seen already the non-zero elements of ( m must be Q, Cj_5 etc, and to have 
discrete but not continuous symmetry at least two of such elements must be non-zero. For 
L = 1 and 2 this is not possible. For L = 3, the only possible pair is (£3, £-2) if j = 3, 
(equivalently (£ 2 , C-3))) an d for L = 4 an additional possibility (£4, C-i) if j = 4. For these 
states to be inert one must have other symmetry operations in H that would impose a rela- 
tion between these non-zero elements. Clearly a horizontal rotation (and hence groups like 
D 5 ) is not sufficient, and hence one is left with the icosahedral group Y. We recall 25] first 
that the icosahedral group Y contains, in additional to the five-fold axis pointing along the 
(12) vertices, three-fold axes perpendicular to the (20) triangular faces and (15) two-fold 
axes pointing out to the (30) edges of the icosahedron. Recall also that in the group Y, 
we have C5C3 = C2, where C3 is a three-fold axis normal to the triangle which contains 
the C5 as an vertex, and C2 is around the edge of the triangle connected to the mentioned 



vertex for the axis 



251 ] . In the isotropy group H, these two- and three-fold rotations C2 



and C3 must appear as symmetry elements, though they can be combined with again phase 
factors and/or time reversal symmetry. By similar argument as given already for C5, C3 
cannot be appear in the form C^e^G (else the state is again time- reversal invariant), thus 
the three-fold rotations can appear only via C^e 2 - 7 ' 27 ^ 3 for some integer j'. Therefore now H 
must contain the element C 5 e l ^ 2 ' K ^C^e^" 2 ' K ^ i = C 2 e^5 + ^"^ 27r = h 2 and hence h\. But this is 
possible only if j = j' = since Cf = 1. Thus we arrive again at the conclusion that H must 
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contain C5 itself, which however is impossible for L < 4. This completes our proof. 

APPENDIX B: RELATION BETWEEN STATES D AND D' 

We here show that the state xyz, with symmetry O(T) x6 (eq (j3J)), is related to 

D' = (\/2, 0, 0, i\/5, 0, 0, V2) (Bl) 

in a new set of orthonormal axes (u,v,w). Since D' has explicitly three-fold symmetry 
along its quantization axis (identified as w) and two-fold axes (up to e ?7r ) perpendicular 
to it, and xyz has a three-fold axis along the cube diagonal, this suggests that we take 
w = (x + y + z)/V3, u = (—x + y)/\/2 along one of the L^e" 1 " axis, and v = w xu completing 
the triad. This transformation between x, y, z and u, v, w is the same as the one used in the 
5 = 2 cyclic state. We have the relations x = ~-^ u ~ ^ v + V = ~j2 U ~~ ~Jq v + "75 w i 
and z = yff + ~^ w } therefore 

xyz = -V2v(3u 2 - v 2 ) + w(2w 2 - 3u 2 - 3v 2 ) . 



On the other hand, we have (in the u, v, w axis) Y^ 3 = =)=y ^-^(u±iv) 3 , Y® = y ~(5w 2 — 3), 
so state D' is proportional to 

®Li[-y/2v(toi 2 - v 2 ) + w(5w 2 - 3)] 

7T 4 

Using u 2 + v 2 + w 2 = 1, we see that state D' is proportional to i times the state xyz. We 
conclude that D' in (u, v, w) axes is the same as zD = (0, 1, 0, 0, 0, —1, 0) in (x, y, z) axes. 




APPENDIX C: SPIN WAVEFUNCTIONS FOR STATES O AND C OF 5 = 4 

We would like to find the spin-wavefunction corresponding to states O [eq ([IB]) ] and C 
[eq (I2ip ], We note that both these expressions contain only even powers of x, y and z, 
they must be just linear combinations of Y^ + l^ -4 , Y^ + Y^ 2 and Y%. With c 4 = yj;, 

we have F 4 4 + F 4 " 4 = c 4 y/§§ ((x 2 + y 2 ) 2 - 8x 2 y 2 ), Y 2 + Y 4 ~ 2 = c 4 ^j(7z 2 - l)(x 2 - y 2 ) and 
= c 4 |(35;z 4 — 30z 2 + 3). For state O, we note further that eq (jTBJ) is symmetric under 
x 2 <-> y 2 , thus Y 2 + Y^ 2 is not involved (as expected since its isotropy group contains the 
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element C4). Noting that Y 4 ° is cylindrically symmetric but Y 4 + Y 4 4 is not, we find easily 
that the expression eq (ll~8l) can be rewritten as 

-> + + 

and hence the spin-wavefunction eq (fl~9l) . 

For the state C, it is more convenient to work with 

vl/' = 7(«/V + e W + x 2 y 2 ) + (ex 2 + e 2 y 2 + z 2 ) . (CI) 

Using the definition of e, we get, 



„2 „.2\ , 2 



'(C2) 

Our task is again simplified by noting that Y 2 + Y^ 2 is the only term that is asymmetric 
under x <-> y, Y^ + Y^ 4 is the only term that is symmetric under x <-> y yet not cylindrically 
symmetric. We can then find easily 

v = ^4 - yf *(n 2 + >7 2 ) - / J(n 4 + n" 4 )] (C3) 



hence the wavefunction is proportional to (— y^j, 0, — y|i, 0, 1, 0, — yfi, 0, — y^j) and 
therefore the expression ( 122|) as shown in text. 
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